Scanning of the e + e — > 7r + 7r cross-section below 1 GeV 
by radiative events with untagged photon 

V.A. Khoze 1 , M.I. Konchatnij 2 , N.P. Merenkov 2 , G. Pancheri 3 , 
L. Trentadue 4 , O.N. Shekhovtzova 2 



1 Department of Physics and Institute for Particle Physics Phenomenology, University of Durham, 
DH1 3LE, UK and Petersburg Nuclear Physics Institute, Gatchina, 188350, Russia 

2 National Science Centre Kharkov Institute of Physics and Technology, 61108 Akademicheskaya 
1, Kharkov, Ukraine 

3 INFN Laboratori Nazionali di Frascati, P.O. Box 13, 00044 Frascati, Italy 

4 Dipartimento di Fisica, Universita' di Parma and INFN, Gruppo Collegato di Parma, 43100 
Parma, Italy 



Abstract 

We discuss an inclusive approach to the measurement of the e + e~ — ► ir + ir~ cross- 
section by the radiative return method without photon tagging. The essential part of 
this approach is the choice of rules for event selection which provide rejection of events 
with 3 (or more) pions and decrease the final-state radiation background. The radiative 
corrections to the initial-state radiation process are computed for DA<1>NE conditions, 
using the quasi-real electron approximation for both, the cross-section and the underlying 
kinematics. The two cases of restricted and unrestricted pion phase space are considered. 
Some numerical calculations illustrate our analytical results. 

PACS: 12.20.-m, 13.40.-f, 13.60.-Hb, 
13.88.+e 



1 Introduction 



The recent high precision measurement of the muon anomalous magnetic moment (g — 2) M Jlj 
has boosted interest in renewed theoretical calculations of this quantity , since any difference 
between the experimental value and the theoretical evaluations based on the Standard Model 
(SM) may open a window into possible new physics ||. While conclusions about posssible 
discrepancy with the SM are premature || , the Brookhaven based experiment is now planning 
a new measurement with three times better accuracy, which may create further challenges to 
the theory. Presently, there are two main sources of theoretical uncertainty in the calculation, 
namely the impact of the light-by- light contribution @, ]7J and the estimate of the error 
from the hadronic vacuum polarization contribution to (g — 2)^. In this paper we address the 
question of this error, for which different groups give different results [§, f|. 
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The problem of the hadronic vacuum polarization contribution is that it cannot be calcu- 
lated analytically because perturbative QCD loses its predictive power at low and intermediate 
energies, where, on the other hand, the effect is the largest. However one can evaluate this 
hadronic term from the data on electron-positron annihilation into hadrons by using a disper- 
sion relation The necessary condition for a theoretical error matching the experimental 
accuracy reached in the (g — 2) M measurement, is the knowledge of the total hadronic cross 
section with better than one per cent accuracy. The recent precision measurements of the total 



hadronic cross-section by the CMD-2 |LT] and BESII [0 collaborations were included in the 
new analysis of Refs. [jTS], O]. While this reduces the error in the hadronic contribution to 
the shift in the running electromagnetic coupling, for the muon (g-2) value it is mandatory to 
perform new measurements of the total cross section at energies below 1.4 GeV (in particular, 
in the e + e~ — > 7r + 7r~- channel) with at least one per cent accuracy. Such accurate measurement 
will then be important not just for the muon anomalous magnetic moment but also for testing 
the effective fine structure constant. 

In the last years, the idea to use radiative events in electron-positron collisions for scanning of 
the total hadronic cross section has become quite attractive. The radiative return approach was 
first discussed long ago, and the lowest-order cross sections for the radiative process of electron- 
positron annihilation into a pair of charged fermions or scalar bosons were calculated ||15|. This 



subject was subsequently studied in several papers (see, for example, [16], [17], 0, [19], P0|)i where 
higher-order radiative corrections were taken into account. 

Due to differences in the systematic uncertainties in the measurement, the radiative return 
approach has several advantages when compared to the conventional energy scan : for example, 
luminosity and beam energy effects are accounted for only once. For variable total hadronic 
energies from the Im^ threshold up to 1.02 GeV, the ideal machine for scanning the total cross 
section, using the radiative return method, appears to be the DA$NE accelerator, operating 
at the $ resonance, together with the KLOE detector [pi], ^2|. DA$NE measurements can 



become quite competitive to the conventional direct cross section scan, and, as mentioned, 
have certain advantages due to the systematics. The radiative return method allows to per- 
form precise measurements of the hadronic cross sections in the p resonance region. The high 
accuracy of scanning is provided by the high resolution measurement of the pion 3-momenta 
(and consequently the invariant mass) with the KLOE drift chamber. Recently the first prelim- 
inary results of such measurements of 7r + 7r _ production cross section below 1 GeV have been 
reported |22[ . 

In our previous paper |20|] the analysis of Initial State Radiation (ISR) effects, which provide 
the basis for the radiative return strategy, has been performed for the realistic conditions of the 
KLOE detector. It was assumed that both the energy of the photon in the calorimeter, and 
the invariant mass of the 7r + 7r~ - system were measured. 

Notice that the KLOE detector allows to register photons only outside two narrow cones 
along the beam directions (the so-called blind zones ). Because of this geometrical restriction, 
most of the ISR events become inaccessible for tagging and cannot be recorded by the photon 
detector. This decreases the statistics and, thus, results in lesser precision. In order to avoid 
this problem and, moreover, to fully exploit the possibility of high precision measurement of 
the two charged pions with the drift chamber [18|, |21|, it was proposed^ to make the photon 
tagging redundant. The idea is to use an Inclusive Event Selection (IES) approach, in which 
only the invariant mass of the final pions is measured, and the ISR photon remains untagged. 

As briefly discussed in [20|, one of the main advantages of the IES strategy is the rise of the 



1 Our attention was first drawn to this idea by G.Venanzoni (see also |20|| ) 
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corresponding cross section caused by the ln(E 2 /m 2 ) enhancement (here E is the beam energy 
and m is the electron mass) due to the possibility to include events with ISR collinear photons 



which, otherwise, belong to the blind zones. As shown in f23fl , the number of such untagged 
photon events exceeds by about a factor three the number of events with the tagged photon (if 
the opening angle of the blind zone equals to 10°). 

Of course, in order to avoid uncertainties in the interpretation of IES approach and to have 
the possibility to describe IES in terms of ISR events, some additional event selection criteria 
should be imposed. The corresponding additional restrictions should make the photon tagging 
redundant, but at the same time should guarantee the suppression of the main background 
caused by events from $ — > 7r + 7r~7r° decay. 

In this paper we present the analytical calculation of the Born cross section of ISR process 

e~(pi) + e + (p 2 ) -> 7 (fc) + tt+(p + ) + 7T(p_) (1) 

and the QED radiative corrections (RC) to it for IES setup, accounting for the additional 
kinematical constraints on the event selection, which can be realized at KLOE. The physics 
motivation for these constraints is discussed in Section 2. 

The Born IES cross section is calculated in Section 3. Note that for a chosen set of selection 
rules, IES cross section at the Born level coincides with the tagged photon events cross section 
as given in Ref. [TjJ, provided that the final pion phase space is unrestricted. But we consider 



also the realistic case when the pion phase space is restricted. 

In Section 4 we discuss the RC to the Born cross section caused by the emission of real and 
virtual photons. At the RC level, the IES cross section differs from the tagged photon result 
because of the contribution of double photon bremsstrahlung. The situation here is similar 
to the case of radiative corrections in DIS with detected lepton (the analogue of the tagged 
photon events) or hadrons (analogous to the IES). In the latter case, the radiative corrections 
factorize while in the former they include by necessity some integrals over hadronic cross section 
that have to be extracted from experimental data. This fact certainly makes the IES approach 
more advantageous. In Section 5 the cancellation of the infrared and collinear parameters, used 
in the calculations of radiative corrections, is demonstrated, and the expression for the total 
photonic contribution to the radiative corrections is given. In Section 6 we discuss also possible 
contribution of the e + e~-pair production into the IES cross section if the e + e~7r + 7r~ final state 
is not rejected by the analysis procedure. Our Conclusion contains a brief summary and the 
discussion of the background processes which may contribute into the IES cross section. 



2 IES selection rules 



As mentioned in the Introduction, the main condition of the IES approach is the precise 
measurement of the di-pion invariant mass in process (1). In addition, restrictions must be 
imposed in order to select final states with only 7r + 7r~ + wy , excluding 7r + 7r~7r°. Finally, we 
have to add some constraints in order to reduce contributions from final state radiation (FSR). 
As an added bonus, such contraints also simplify the theoretical calculation of the radiative 
corrections. 

Rejection of the 3-pion final state in process (1) can be done selecting events with an 
appropriately small difference between the lost (undetected) energy and the modulus of the 
lost 3-momentum in process (1). In terms of the measured pion 3-momenta, this restriction 
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reads IffSl EDI EI 



2E-E + -E_- |P $ -p+-p_| , 77 «1, (2) 

where E is the beam energy, E± = yp± — m 2 is the energy of 7r ± , and is the pion 
mass. P<j> is the total initial-state 3-momentum which, at DA$NE, is non-zero, due to a 
small acollinearity in the beams, |P$| = 12.5MeV. If the chosen parameter 77 is small enough 
(< ( m ir/E) 2 ), constraint (2) allows to avoid the undetected tt° and to retain only the undetected 
wy system. Inequality (2) can be rewritten in terms of the total energy Q and modulus of the 
total 3-momentum |K| of all photons in the reaction e + + — ► 7r + + 7r~ + 77,7 as 

^ - |K| < 77^. (3) 



The optimal value 77 = 0.02 decreases also the FSR background |L8[ . 

The next constraint selects such events, where at n — 1 the undetected photon is collinear 
with the emitting electron (or positron). The collinear events considered here are those in which 
the photon belongs to a narrow cone with opening angle 26o (9o <C 1) along the electron beam 
direction (the blind zone for the KLOE detector). This constraint reads 

Kpi > \K\Eco, c = cos# , (4) 

where pi is the 3-momentum of the electron and 80 can be chosen to be 5° 4- 6°. Due to this 
constraint, the collinear photon radiated by the initial electron contributes to the observed IES 
cross section and induces a ln^^/m 2 ) enhancement, which, at DA$NE, makes the IES cross 
section a few times larger than in the tagged photon case. Moreover, the collinear constraint 



provides the possibility to apply the well known quasi-real electron (QRE) method [24j to 
calculate radiative corrections. Even at the Born level, the difference between the exact result 
and the corresponding QRE approximation is negligible (see for details Section 3). Selection 
rules (2), (3) and (4) imply precise measurements of the pion 3-momentum that can be provided 
by the KLOE drift chamber. 

Note also that in the Born approximation (n = 1) inequality (3) is always satisfied, and, 
therefore, the selection rules (3) and (4) imply non-trivial consequences only through the con- 
tribution to the radiative corrections from two hard photon emission. 

Because of the existence of the blind zones, the KLOE detector cannot provide the detection 
of the final 7r + and tt~ inside the full phase space, picking out events with pion polar angles in 
the region 

e m < o± < 71 - e m . (5) 

In principle, 8 m can be taken to be about 10°, but, as shown by the Monte Carlo calculations 



T% jlBf, the choice of 6 m influences also the value of the FRS background. The optimal value 



of 9 m for DA$NE conditions is 20°. 

Usually, the restricted pion phase space can be taken into account by introducing an accep- 
tance factor A(9 m ). The calculation of this factor is very simple for a nonradiative process, but 
for the ISR process, and the RC to its cross section, it is non-trivial. The analytical form of 
A(6 m ) can be derived in the framework of QRE approximation. To illustrate the problem, in 
the following we perform calculations for both unrestricted and restricted pion phase space. 
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3 Born approximation 



To lowest order in a, the differential cross section of process (1) can be written in terms of 
the leptonic and hadronic tensors as (see Ref. 

B M 2 a 2 a d 3 k d 3 p + d 3 p_ 

*> = ^ L UP^m^ 2 — 16%2E+E J (<l-P + -P-) , (6) 

where us is the energy of photon, and the hadronic tensor is expressed via the pion electromag- 
netic form factor F n (q 2 ) as follows 

H^ u = -M F Al 2 )\ 2 P-nP-» , P-h = P-h-Iqh, q = Pi+P2~k = p + +p_ . 

The pion electromagnetic form factor defines the total cross section u{q 2 ) of the process 
e + + e~ — > 7r + + 7r~, by means of relation 



|F. W = C = (l - • (7) 



In the case of unrestricted (full) pion phase space, the integration of the hadronic tensor can 
be performed in invariant form [U 



167T 1 



f d 3 p + d 3 p„ q 4 2 _ q^q u 

J ^ E £L ( Q ~ P+ ~ P ^ = 8^ a(g 9flu = 9veu ~ f ' ® 



whereas for the restricted case we have first to contract the tensors in Eq. (6) and then to 
integrate the result over the pion phase space. 

The leptonic tensor on the right-hand side of Eq. (6) for the case of collinear ISR along the 
electron beam direction is well known [15|, 25 1 

T ~ , r \ \{q 2 -U) 2 + {q 2 -t 2 f 2m 2 q 2 i Aq 2 _ „ / Aq 2 8m 2 \ „ 

k) = — "2— g„ v + —Pi^PXu + [ — " 72- )P2uP2u, 

(9) 

where 

Pi,2|i = Pi,2» ~ ' *i = ~ 2k Pi, h = ~2kp 2 , s = 2p l p 2 , q 2 = s + h+ t 2 , 

and we neglect terms of the order \m 2 /t 2 \ which are always below the required accuracy.^ The 
contraction of the leptonic tensor with that is necessary to use in the case of unrestricted 
pion phase space reads 



r7 / ,,~ 9F 771 (q 2 -ti) 2 + (q 2 -h) 2 2q 2 m 2 

\Pi,P2, kjg^u = IF , F = — tj • ( 10 ) 

iii 2 t l 



In the case of restricted phase space we have to use the following relation 

LU Pl , P2 , k)H. v = Aq 2 \F n (q 2 )\ 2 R , (11) 



2 Note that formula (9) corresponds to the radiation along the electron beam direction (see also inequality 
(4)). To account for the radiation along the positron we have to multiply all results for the cross sections by a 
factor of 2 (see the end of Section 6). 
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m 



2 



Q 2 (Xi + X2) ~x\ -X2 Xi X2 , 2m 2 x 2 /X2 



""7T 7^ 1 V VA1 1 A,zy A.1 A2 A.1 AZ . — AZ / A^ -1 \ 

<T M2 h t 2 t{ \q 2 J 

Therefore, the differential cross section of the process (1) in the Born approximation can be 
written as 

drf = atf^F^ , (12) 



da 



R 



3a(q 2 )R a d 3 k \p_\dE_dc_ d<p_ 
< ~E~+ 2T 



■5(2E - uo - E + - EJ) , 



(13) 



where c_ = cos 9^ and 0_, (p_ are polar and azimuthal angles of the negative pion respectively 
(we take Z axis along the pi direction). 

It seems at first sight that one can perform the trivial integration with respect to the 
azimuthal angle ip- on the right-hand side of Eq. (13) because the quantity R does not contain 
any ^.-depending term. But in the general case the pion energies E + and E- depend on </?_. 
Moreover, the upper limit of integration over c_ for the ISR events with collinear photon along 
the electron beam direction is smaller than c m = cos9 m and depends on </?_ as well, i.e. 



< C r 



(14) 



where c max must be determined from the 3-momentum conservation, provided that c + = 
cos# + = — c m . Therefore, in the case of the restricted pion phase space, it is necessary to 
first integrate over c_ and then over </?_. 

Let us now perform the integration with respect to the photon angular phase space on the 
right-hand side of Eq. (12). It is convenient to choose X axis along the P$ direction. In the 
laboratory frame 



p2 



pi = (£,0,0,|pi|) , p 2 = (f?,|P*|,0,-p a? ) , P2z = E(l-^) 



(15) 



and 



ti = -2u(E- |pi|cos^) , t 2 = -2uE[l+ (l - cos 9] + 2cu\P^\ sin cos 0] , 

2E 



s = 4E7 2 -P|, q 2 = 4E 2 -4Euj-Pl(l-^cos9) + 2cj|P $ | sin# cos0 , 



(16) 



where 9 and <fi are polar and azimuthal angles of the photon radiated in the initial state. 
Keeping terms of the order Y\jE 2 and 9%, we can use the following list of integrals 



r do n9 2 r do 

J ~ 2 ~ 4oj E ' J ~i 



00 qE 



4E 2 J 



2-, 

12 



r dOm 2 ti 

' J = ^ 



p2 . 
*. 1 



2E 2 J ' 



/ 



dOt 2 



= 47T 



p2 

4E 2 , 



x 6»g-| /-dO 



7T 



4^ 2 £ 2 



AE^) Lo + J 



(17) 



4£ 2 



CJ 



4£ 



p2 E 2 9 2 
— - , L = In — , dO = d cos 
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When evaluating these integrals we systematically neglected small terms of order ~P\9 2 /E 2 . 
Within such approximation we can use the substitution 

|P$|sin0cos0 -> 0, P|cos0-+P|, (18) 

in the expressions for the invariants t\ and £2 in Eqs. (17). With the same accuracy we can 
write 

ujduj 



^Ml + 4). (19) 
4E \ 2E 2 J v ; 



Combining (12), (17) and (19), we arrive at the distribution over the pion squared invariant 
mass q 2 , for unrestricted pion phase space 



dap _ a(q 2 ) a \( q 4 , ^ , u °\(i \ P * V ^ , °o ( ^ ] Q 2 2uJ c 



dq 2 AE 2 2TiV\SojQE IE 2 E J \ AE 2 J u 2lu E 6 \8u E 2E 2 E 



(20) 



To guarantee only one per cent accuracy one can neglect terms proportional to P|/4E 2 and 
#q/6 in (20) because their contribution into the IES cross section is of the relative order 10 -4 . 



Such procedure leads to the well known result corresponding to the QRE approximation [23 



dcr F _a(q 2 )a _ 2z q 2 

W ~ l£^ P(z ' Lo) ' P{z > Lo) ~ T^7 L ° " — ' z ~ ae 2 ■ (21) 

Thus, the QRE approximation appears to be sufficient for a description of the IES cross 
section and corresponds to a one per cent precision even at the Born level. 

Consider now the case of the restricted pion phase space. First, the 5-function on the 
right-hand side of Eq. (13) has to be used to perform the integration with respect to E_. Then 
within chosen accuracy we obtain 

|P_|,/ B ~^ E -<*-E + -E-.) = ,„ , I ' P :| 2 ^7— 3T^T , (22) 



E. 
where 



( P 2 \ 

c 7 _ = cos0c_ + sin$,s„ cos ((f) — (pj) , s_ = sin^. , u> — uj [1 + ^^rj • 

Now let us express E_ in terms of the photon energy uj and angles of the photon and negative 
pion, using the energy-momentum conservation. The result can be written in the following 
form 



*- = AB - C %;^- C2) , * = B = 4E( B - W )-|P.I«, (23) 

C = UJC 1 _ — |P$|s_ COS<y?_ . 

To find the energy of the positive pion it is necessary to use the relation 2E — oj = E + + E_. 
Having expressions for the pion energies, we can apply conservation of the Z-component of 
3-momentum at c + = —c m 

— |p+|c m + |p-|c_ + ucos9 = (24) 
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to derive the upper limit c max of the variable c_. In Eq. (24) both |p_| and |p + | are functions 
of uj, c_, cos 8 and cos (0 — </?_). Therefore, to reach the same accuracy as in Eq. (20) for the 
case of the restricted pion phase space, the differential distribution over the pion invariant mass 
squared has to be taken in the form 



B i 2-7T 2w C max 

da R _ 12a(g 2 ) a ow P| \ f [ f d<p. R\p^\ 2 

dq 2 ~ AEX 2nE{ 1 + 2Ei) J dC0S& J 2nJ 2tt J & C ~ A|p_| + CE. ' {2b) 





2 




P-l 





CO 



where Co = cos#o and the quantities A, C are defined in (23). The analytical integration on the 
right-hand side of Eq. (25) is not available and the task of integration can be left for numerical 
calculation. 

The result is very much simplified if one neglects terms of the order |P||/i? 2 and use the 
QRE approximation, assuming cos 9 = 1 in expressions (22), (23) and (24). This leads to the 
IES cross section 

da R cr(q 2 ) a „. , . . . . 

5? = ^S P ^^-^' (26) 

where z is defined in (21) and 

A(7 v_ 12 °T Hr n z[(l + z)K - (1 - z)c_f TJ _X^_xL_rf L _ (97] 
l ' Cmj "C J K[(l + zf - (1 - zfcV? ' U ~ 4E* 16£ 4 4zE* {<} 

represents the acceptance factor (see the end of the previous Section) as it follows from the 
comparison between the IES cross sections (21) and (26). To write down A(z, c m ) the following 
notations were used: 



21 



4E 2 (1 + zf - (1 - z) 2 cl ' 4E 2 ' 



/ n_ (1 + z)g _ zc m \(\ +z)K(cm) + (1 ~ z)cm] 1-z 

° max[Z ' ° m) ~ y/(z - (1 - z)gf - (1 + z)W ' 9 = ~ (1 + z? - (1 - zfcl ~ 2 • 

(28) 

At fixed values of c m , the quantity c max depends on the squared pion invariant mass q 2 = 
AE 2 z. If z is small enough (this situation corresponds to the radiation of a very hard collinear 
photon with the energy fraction (1 — z) by the electron) c max , as formally defined by (28), can 
approach — c m and became even smaller. Because selection rule (5) forbids any values of c max 
smaller than c m , it is necessary to substitute the upper limit of integration on the right-hand 
side of Eq. (27) by max[c max , —c m ]. If not, the formal calculation in (27) leads to negative 
values for the acceptance factor at small z, as one can see from from Fig. 1 for c m = cos 20° 
and z « 0.1, while really it equals to zero at such ^-values. Note that in the framework of the 
QRE approximation one can use also energy conservation for events with 6 + = n — 9 m 

E(l + z) = E_(c max ) + E + (c + = -Cm) 

in order to obtain the analytical form of c max on the Born level. 
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If c m = 1, 3-momentum conservation in process (1) requires c max (z, 1) = 1 as well. In this 



case 



K(Cr, 



A5 2 

z 



[l + z)K(Cm)-l + Z 

4 



and it is easy to see that c max (z, c m ), as given by (28), satisfies this requirement. 

In principle, the acceptance factor A(z, c m ) may be computed analytically by means of 
Euler's substitution on the right-hand side of Eq.(27) 



2t 



K 



5(1 - z)(t 2 + a 2 ) 
2zt ' 



dc. 



(t 2 + a 2 )dt 
2t 2 ' 



S 2 (l + z) 



5 2 (l-z) 2 



By definition, A(z, 1) = 1 and for c m < 1 always A(z,m) < 1. 

We leave the task of the analytical integration of the acceptance factor in general case aside 
and first only note that the limit 5^0 (that is, of course, not the case for DA$NE) can be 
used to control our calculations. In this limiting case 



A(z,K = 1,1) 



12z 2 



dx{l — x 2 



1 + z + fl - z)x\ 



1 . 



The acceptance factor as a function of the pion squared invariant mass is shown in Fig. 1 
for 9 m = 10° and 9 m = 20°. We see that the acceptance factor A(z, c m ) is close to unity in a 
wide z-range, but decreases very rapidly with the pion invariant mass. 



A(z,cj 
1.00 - 



c =cos(10 ) 



0.0 0.2 



A(z,c 



c =cos(20 ) 



Figure 1: Acceptance factor defined by Eq. (27) for different values of c m 

Finally, let us note that two representations (21) and (26) for the IES cross section can be 
obtained also by inserting the QRE form of the leptonic tensor 

a d?k a dz 2 ,1_ 2z 2 _ _ , 

4:7i z uo p 2tt z 2 q z 

into Eq. (6). 
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4 Radiative corrections 



If events with e + e~7r + 7r~ final state are rejected, only photonic RC have to be taken into 
account. These corrections include contributions due to virtual and real soft and hard photon 
emission. To calculate them we use the QRE approximation from the very beginning. 



4.1 Soft and virtual corrections 

The soft and virtual corrections are the same for both, unrestricted and restricted pion phase 
space, and the corresponding contribution can be found by the simple substitution 

£-P(z, Lo) - (Pj 2 C S+V , C s+V = P P(z, Lo) + D(L S , L , z) + N(z) (29) 

Z7T \llxJ 

in the right-hand sides of Eqs. (21) and (26). As a result, we have 

, S+V , o\ g , S+V j s+v 

d(TF -°b)(lL\ 2 C 8 +v , ^ = ^(^1. (30) 



dq 2 AE 2 \2nJ ' dq 2 dq 2 

All the logarithmically enhanced contributions to C s+V are contained in the first two 
terms in Eq.(29) and were first found in Ref. |16[-0 The third term in (29) describes the 
non-logarithmic contributions [2"B||. Finally, we obtain 

2tt 2 9 AE 2 
p = A(L S -l)ln A + 3(L s + \nz) + — -- ,L, = hi— , (31) 

3 2 m A 

1 + z 2 7r 2 1 + 2z — z 2 Az In z 
D(L s ,L ,z) = L [(L -2L S - In z)\nz + — -2Li 2 (z)] + — —L + L s , 

, 1 8,2 In z r o ln 2 z-, 7r 2 . 5 N lt „ 5 w . . . 

N(z) = — + 2z In 2 (1-z ) + - + — (Az + 6-- ) + (4z-6 + - )Li 2 (z) , 

1 — z 1 — z L 1 — z 1 6 1 — 2 /v 1 — z' 

where AE (A <C 1) is the maximum energy of a soft photon, and 

X 

Li 2 (x) = - [ — ln(l - y) . 

J y 



4.2 Two hard photon emission along the electron and positron di- 
rection 

Concerning the contribution from additional hard photon emission, we divide it into three 
pieces. The first one is responsible for the radiation of an additional photon with energy 
uj 2 along the positron beam direction (provided that a collinear photon with energy u)\ is 
emitted along the electron beam direction). To calculate it we introduce the angular auxiliary 
parameter fljcl and use the QRE approximation to describe the radiation of both photons. 
The calculations are the same for restricted and unrestricted pion phase space and do not affect 

3 There is a misprint in the expression for p given in Ref. pq| . The correct expression includes an additional 
7r 2 term. 
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the acceptance factor (as in the case of virtual and soft corrections). Using the subscript F to 
indicate the cross-section for the unrestricted phase space, the result reads 

P(xL)P(vL') dy L>-\n E2d '° x-l- Ul v - 1 - ^ 
~W~^&L ^ P{x,L )P{y,L )j, L -ln— ^, x-l -, y - 1 -, 

(32) 

where q 2 = 4E 2 xy is fixed. The maximum value of the photon energy uj 2 can be obtained 
from restriction (3), taking into account that, for the events under consideration, the QRE 
kinematics (cos6*i = 1, cos^ = —1, 9i j2 are polar angles of photons) gives 

Q = u 1 + uj 2 , |K| = ui - u 2 , -> lo 2 < — , y Q = l-~. (33) 

Because of the smallness of the parameter 77, only terms singular at y = 1 contribute to the 
integral on the right-hand side of Eq. (32). Moreover, we can substitute x with z, neglecting 
terms of order rj as compared with unity (because q 2 = 4E 2 z = 4E 2 xy), and perform an 
elementary integration over y, obtaining 

4.3 Two hard photons emitted along the electron directions 

The second contribution to the radiative corrections caused by an additional hard photon emis- 
sion corresponds to the radiation of two hard collinear photons (each of them with energy larger 
than AE) by the electron, provided both are emitted within the narrow cone of opening angle 
26q along the electron beam direction. This contribution also does not affect the acceptance 
factor, and the result can be written with the same accuracy as (30) and (32) 

da" F _ a(q 2 ) f a\ 2 h da" R _ da" F 

C H = B x (z, A)L 2 Q + B 2 (z, A)L + B 3 (z, A) . 
Functions B\(z,A) and B 2 (z, A) were first calculated in Ref. |2j 



B 1 (z, A) = l -P 2e (z) + (]nz - ? - 2 In A) , 

/an / x (3 + z 2 ) In 2 z 2(1 + z) 2 , 1-z 
fl,(»,A)=3(l-,) + i 2{l ! :) -- L—Lln— , (36) 

where P 2 g(z) is the ^-term of the second order electron structure function (see, for example, 

'^ + z 2 /, s , 3\ 1 



(21n(l-z) -lnz + ~) + -(1 + z) In z - 1 + z 



1-z 

Function B 3 (z, A) was calculated in Ref. [EBJ and reads 



n . A . Az 1-z Az tt 2 Az ( 17 28 8 

B 3 (z,A) = n — ; h V —, r ; nz+ (37) 

3V ' ; 1-z A 3(1 -z) 6 3 V 3 3(1- z) 3(1 -z) 2 J 
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J 



3-2Az + 5Az 2 -A8z 3 + 7z\ 2 . 7*\ ,-•/-. s T 

+ In 2 z+(l- —)Li 2 (l -z) + J, 

1-2 1 

z 2 + (1 - x) 4 f dt _ , . z + x , xA - 3,2 

— ry- rj- / — L(t,x,z +7771 V Ll x ' z +^71 ^L 2 {x,z) 

xX{l — xy J t 2(1 — x) 2(1 — x) z 





where the logarithmic functions L, Li and L 2 are defined as follows 

(1 - x) v/F(t, x, i) + tx(* - xX) + A(l - x) 2 

L(t, x, 2) = m 1 ? r-T — , A = 1 — x — z, 

2(l + £)A(l-s)* 

r , , (x + z))JF(l,x,z) + X(z - xA) + x(x + z) 2 
U (x, z) = In — , 

Tl , . (1 -x)a/F(1,x,;z) + x(z-xA) + A(l - x) 2 
L2(X ' Z) = ln 2A(T^ ' 

F(t, x, z) = A 2 (l - x) 2 + 2txX(z - xA) + t 2 x 2 (x + zf . 

The integral over x in the expression for J diverges when x — > and x — > 1 — z. But the 
structure of the integral is such that these divergences compensate each other, and this can be 
seen by taking into account that: 

i) integral overt converges, 

z 2 + (l-x)\ u 1 + z 2 t(l-z) 
ii) lim — — - — L{t,x, z) — — -In- 



x^o xX(l — x) 2 ' ' x(l — z) xz 

z 2 + (l-x) 4 r/ . 1 + z 2 , t(l-z) 

The quantity J can then be computed numerically with appropriate precision. Note also 
that the integral over t on the right-hand side of Eq. (37) can be computed analytically (see 
also i.e. 

i 

A 



f dt _ . . 1 9 . . _ . (2xz\ _ . / 2x 2 A\ _ . / 
/ -y L(f, x, 2;) = -L 2 (x, z) + Lz 2 (— j + L% 2 y —j + Lz 2 (- 



o 



/l = A(l - x) 2 + x(z - xA) + (1 - x)y/F(l,x,z) . 

4.4 Emission of one hard collinear and one hard wide angle photon 
with unrestricted pion phase space 

The third, least trivial, contribution into the RC caused by two hard photon emission is con- 
nected with events when one photon with energy ui is collinear and the other (with energy 
u 2 ) is radiated at angles between 7r — 9' and 6> . For such events, constraints (3) and (4) are 
somewhat tangled, and one needs to choose convenient variables to disentangle them and to 
determine the photon phase space. In spite of the obvious fact that these constraints concern 
only the photons, their respective contribution to the RC affects the acceptance factor given by 
Eq. (27). The physical reason for such a behaviour is that now (in contrast with previous cases) 
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the 3-momenta of photons and pions do not lie in the same plane even if the QRE kinematics 
(ki = (1 — x)pi) is applied. 

Nevertheless, the QRE approach for the description of collinear photons allows to simplify 
the form of the cross section and to disentangle all the kinematical restrictions. In accordance 
with this approach, the starting point for our calculation of the differential cross section, suitable 
for the unrestricted pion phase space case, is the following 

H a(q 2 ) a a dxd 3 k 2 , . 

da * F = ^^ P ( X ' Lo)L ^ XPuP2 ' k2)9 ^4^^J^ ' (39) 

= 2lTUJ 2 duJ 2 dc 2 , c 2 = COS #2 , x — 1 - ^ , 
UJ 2 E 

where 9 2 is the polar angle of the non-collinear photon. Since our aim is to derive the differential 
distribution in the squared pion invariant mass q 2 , it is convenient to use the relation between 
q 2 and c 2 to avoid the integration over c 2 on the right-hand side of Eq. (39). In addition, it is 
convenient to introduce the total photon energy fl instead of u 2 

da 2 

q 2 =AE(E-n) +2ccw 2 (l-c 2 ), u 2 = ^ - u u dc 2 -> — — , du 2 = dVL. (40) 

1uj\UJ 2 

In this case in the lepton tensor IS 1 the electron mass can be neglected. Thus, the differential 
cross-section has the following form 

# = 3£(f )><*. *»> [— - ^i- + -) + — + — l -#^} • <«) 

dq z AE Z \2tt/ I Vxu\U 2 \xu\ u 2 J u 2 xui J u>i{E — ui) > 

4E 2 Q Z 4E 
ui = -2k 2 pi = -, u 2 = -2k 2 p 2 = - tvi) - EQ Z ], Q z = Q - E{1 - z) . 

It is useful to rewrite the expression in curly brackets in a form which is convenient for the 
integration over lo\ and Q 

/ T zE 2 L [2z - (1 + z)L }E 
\ L ° " (E-ujJ 2 + E^ + 

[2(1 + z 2 )L -Az- (1 -zf\E 2 - (l-z)(n-2uj l )E^duj 1 dn 
wi(ft-cji) — EVL Z J Ett z ' 

where we neglect terms which do not contain in the denominator the small quantity Q z , which 
is of order r]E, as one can see below from the expressions for VL min and fl ma x- 

Our task now is to define the integration region on the right-hand side of Eq. (41), which 
is determined by restrictions (3) and (4) for the event selection, as well as by the inequalities 

-eg < c 2 < c , EA < ooi < Q - EA, c' = cos6» , (42) 

limiting the possible angles of the non-collinear photon and energies of the collinear one. Re- 
striction (3) defines the maximum value of Q, whereas restriction (4) defines the minimum value 
of u>i at fixed Q 

Q max = E(l - z) (1 + h , oo mm = 2 ™ z , |K| = y/W-4ESl z . (43) 

Z i l — is. Cq 
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To obtain the minimum value of Q we use the first relation in (40) at the minimum possible 
value of u 2 = AE and c 2 = Co 

{l mm A(l - c ) = 2E[Q mm - E{1 - z)\ 

which leads to 

n min = E(l-z)(l+ Ail ~ Co) ) . (44) 
From the condition c 2 > — c n it follows also that 



W - <Wl <o; + , c ± = -[l±^/l--^(l + ^)] . (45) 
Finally, inequality c 2 < c reads: if the values of fl's are such that 



n<n c , n c = E(i-z)(i+ il - Co l il - z) ) 



8 

then 



^ r-i / 8E£l z -, . . 

Wl > ^ + or Wl < w _, w± = -[1 ± ^1 - q2(1 _ Cq) ]- (46) 

To derive the integration region we have to combine consistently all constraints (42)-(46) 
on io\ and Q, and such combination leads to 

u) min < iO\ < and uj + < cu 1 < Q — AE, Q m i n < Q < , = E(l — z)(l + A) , 

^min < u)\ < uj- and u + < u>i < u + , f^A < ^ < fi c , (47) 

The integration region defined by the inequalities (47) is shown in Fig. 2. 

The integration with respect to uj\ and Q on the right-hand side Eq. (40) over the region 
(47) can be performed analytically, and the list of necessary integrals is 

i^^ 1 -^ 2 - 11 - 1 ' ^d-co)(l-,)' 

/ fJ^^n = -In^ln^ + Li 2 (l - z) - Li 2 (-£z) + Li 2 (-£) - Li 2 (-^—^), 
J {h, — u>i)\l z z 

f EduidVl l r . . . _ l + , inS 

J (g-^ -HH^'-d-^— ], (48) 

f (Q - 2uj 1 )du 1 dtt 7T 2 1 2 l-c n (1 - c )(l - c' ) 

7 Mn-.O-^in^ -e"^ — + ln 2A ln I 2Ll ^ 

[ EduidVL + 2 1 ~ c o , V , (l-c )(l-c ) 2 , 

7 l-^ L 2 2 2 2A 4 
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AE 



Figure 2: The integration region with respect to u)\ and Q, as given by inequalities (47) 

Using these integrals we can write the corresponding contribution to the IES cross-section 
from the events presently considered, as follows 



(49) 



G 1 = n -+ 1 -^ 1 -^ 



In^ln^^'^-ln^ 
2A 4 s 



G2 = (1 _ 2)[ _ 2 + hl (l±i|i±i£) ] + 

+ (1 + z) [(1 + InO In ^ + IA 2 (-£z) + L% 2 {- 1 -^-) - Li 2 {-£) - Li 2 (l - z 



Go 



1 - *)(y - ln 2 + 2Li 2 (-0 + 2z[Lz 2 (l 



«)-K a (-^)- J Li 2 (-^)]. 



We would like to emphasize that only the part of the IES cross section, defined by Eq. (49), 
has a non-trivial dependence on the physical parameters 6q and rj, which determine the main 
requirements for the event selection. 



4.5 Emission of one hard collinear and one hard wide angle photon 
with restricted pion phase space 

Consider now the situation with restricted pion phase space. Unfortunately, the calculations in 
this case are not so simple and cannot be performed analytically. Nevertheless, the dependence 
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on the unphysical auxiliary parameters A and 8' , which have to vanish in final result for total 
RC, can be extracted. 

Our starting point is the following representation for the differential cross section corre- 
sponding to the QRE approximation (by analogy with (39)) 

h 12a(g 2 ) a 2 dx _ _ d 3 k 2 dp_ 





P 


2 


(2E-n)\p- 


+ E^ujxC- + uj 1 )c 2 _] 



H = j— -— : — ttt : 7 , C 2 _ = C 2 C- + S 2 S- COS (if 2 - </?_) , s 2 = sin# 2 , 



where ip 2 is the azimuthal angle of the non-collinear photon. 

The experience of previous calculations in the case of unrestricted pion phase space suggests 
that, in order to express the pion energy via angles and photon energies, we can use Q = E(l—z), 
neglecting only small terms of order r\ as compared with unity. Therefore, we have 



E - = 2EZ jt_ ' R = \J l -jM 2 -B 2 ), A = l + z, B = (l-x)c + (x-z)c 2 . 



(51) 



and with the same accuracy 



_ 2z{AK - B) 2 

~ K(A 2 - B 2 ) 2 ' + = ( ' ~ ' ^ ' 

In the limiting case when the energy of the non-collinear photon approaches zero (x — z), 
H becomes twice the corresponding value entering, under integral sign, into the expression for 
the acceptance factor (27) 



x—*z 



The upper limit of the c_ variation has to be determined from the 3-momentum conserva- 
tion, provided that 9 + = n — 9 m , that is 

-|p+|c m + |p_|c_ +E[l-x + (x- z)c 2 ] = , (54) 

where one has to use expressions for E_ and E + to find |p + | and |p_| as given by (51) and 
(52). 

Contracting the indices on the right-hand side of Eq. (50) and using relations (40) we arrive 
at the distribution over the squared di-pion invariant mass 

da" R _ I2a{q 2 ) / a \ 2 p( T t ) duJ ^ dQ TH d{y 2 - , ,,,, 

~d? " ~W P(X ' Lo) (E - Vi "2 2^ dC "' ~E ' (55) 



5 



2 

T = — — 

z 



(q 2 -x Ul ) 2 + (q 2 -u 2 ) 2 i Xi(q -u 2 ) XiW ~ ar«i) x 2 x\ + xl 



XU\U 2 



+ '- + 



U\U 2 XU\U 2 XU\U 2 

To advance further, one must integrate first over c_ because the upper limit c max depends 
on c 2 (as it follows from (54)), and in its turn c 2 is a function of q 2 and of the variables lo\ and Q. 
In the general case we cannot integrate analytically, even to write the analytical expression for 
c max is a problem. But it is necessary to prepare the expressions which can then be integrated 
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numerically. Therefore, the dependence on the unphysical parameters A and 8' has to be 
extracted. 

For this goal, note first that AE is the minimum possible energy for the non-collinear 
photon. Therefore, in order to extract the A-dependence it is sufficient to investigate the limit 
uj 2 — > 0. In this limiting case x — z, and we can use Eq. (53) for the term H. Moreover, 
the 3-momentum conservation (54) becomes the same as in the Born approximation, and its 
solution c max [uj2 — > 0) coincides with the expression given in (28). 

Next, we select terms in the expression for T that are singular in this soft limit, because 
only such terms lead to the A-dependence via In A. From the list of integrals (48) and the 
expression for u\, u 2 (see (41)), it is easy to understand that only the terms containing the 
product u\u 2 in the denominator can induce such a dependence. Taking into account also that, 
in this limiting case, X2 = q 2 — z\i, the necessary terms can be written in the following form 

Tg = 2^16^ 2 = 2E*z(l-z)*U 

U\U 2 £l z [(jJi(£l — UOi) — EVL Z ] 

(for the definition of U see Eq. (27)). To extract the In A-dependence we apply the standard 
subtraction procedure and rewrite the cross section (55) as the sum of its hard and soft parts 



dcr 3R dcr 3R , dcr 3R 



dq 2 dq 2 dq 

where 

da^R _ 12a(q 2 ) / a \ 2 duj^QT.H, 



(57) 



fe) P{Z > Lo) 2E 2 z(l-z) dC -> (K 



dq 2 4E 2 ( \2ii) v ' u; 2E 2 z(l - z) 

and the upper limit of c_ variation in Eq. (58) is defined by Eq. (28). The hard part of the cross 
section on the right-hand side of Eq. (57) is not singular at cj 2 — ■> 0, whereas the integration of 
the soft part over the region (47) induces all A-dependence. Using the corresponding formula 
in the list of integrals (48) we present the soft part of cross section (55) in the form 



da 3R _ a (q 2 



^yP(z,L )G 1 A(z,c m ) , (59) 



dq 2 AE 2 

where G 1 is defined in Eq. (49). 

It is worth noting that this soft part absorbs also all the dependence on the angular auxiliary 
parameter 9' . That is the reason why the hard part of the cross section depends on the physical 
parameters only and can be computed numerically. The soft part (58) of the cross section 
has to be added to the other contributions into the RC to eliminate the dependence on the 
unphysical parameters in the total RC. 



5 Total radiative correction 



The total contribution from radiative corrections to the Born cross section in the case of 
unrestricted pion phase space is represented by the sum 

■j RC , S+V , H , H , H 

aa F ao F da 1F da 2F ^ da 3F 
dq 2 dq 2 dq 2 dq 2 dq 2 ' 
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whereas for the case of restricted phase space we have 



, RC , S+V j H , H , H s , H h 

da R = da R da 1R da 2R da 3R da 3R 
dq 2 dq 2 dq 2 dq 2 dq 2 dq 2 

The two auxiliary parameters, the infrared cut-off A and the collinearity angle 9' Q , enter 
into the individual terms on the right-hand sides of Eq. (60) and (61) in the same combination 

2P(z, Lo) In A [2{L S - 1) - (L' - 1) - (L - 1) + In ^] , (62) 

(here and below we use the expansion of Co and c' ). According to the definition of the large 
logarithms L , L s and L' (see Eqs. (17), (31) and (32)), the expression in square brackets in 
Eq. (62) equals to zero, and, therefore, the total RC depends only on physical parameters and 
can be written, in the case of unrestricted phase space, as 

da R p a(q 2 ) a rc rc a F + L Fi + F 2 

~W = 1^ P(2 ' Lo)6f ' ** = 27 P(z,L ) ' (63) 

1 3 n A 3 n n ^ir 2 Q 

F = -LlP 2e (z)+P(z, L ) [L s (-+2 In |) +ln ^ (-+2 In J-) -2 In 2 £-2 In |+3 lnz+ — - -} , 

3-8z + z 2 2(1 + z) 2 , M , r Az M Wi . I,, 2 l + z 2 tt 2 

-2U 2 (z))+(l-z) In (1 + a . ( 2 1 + ^ ) +(l+z)[£»2(-^)+£» 2 (-— )-£» 2 (-0-£» 2 (l-^)], 

s ^ 

F 2 = (1 - In 2 £ - + — ln(l - z) + [-2* In £ + _ ] In z + 2z ln 2 (l - *)+ 

3 - I2z + 30z 2 - 36z 3 + 7z 4 , , tt 2 , Uz 5 , r , w 5 , 

+ 111 1 + T< 4 + — " — ' + L '^ )i4z " 6 + — ) + 

+(l-i)L» 2 (l-*) + 2Li 2 (-£) -22[Li 2 (-^) + Li 2 (-i^)] + J . 

o Z 

As mentioned in the Introduction, the IES Born cross section and its RC factorize into the 
low energy pion pair production cross-section, u(q 2 ), and a term of a pure electrodynamical 
origin. This latter term depends on the measured pion invariant mass q 2 and on the physical 
parameters rj and 9 which define the rules for the IES. 

The z-dependence of the total radiative correction to the Born cross section (21), 5 F , is 
shown in Fig. 3. Note that the contribution of the non-logarithmic terms to 5 F equals para- 
metrically to a/(27rL ), which is of order 10~ 4 , as the relative contribution of terms proportional 
to P|/4E 2 and 9^/6 in the Born cross section (see Eq. (20)). That is why the exact calculation 
of the Born IES cross-section should be complemented with the radiative corrections calculated 
with the inclusion of non-logarithmic contributions. 

Similarly to (63), we can write the complete radiatively corrected cross-section for the case 
of restricted pion phase space in the form 



da"p cr(q 2 ) a r , . , . c rc da» R 
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Figure 3: The full first order radiative correction to the Born cross section (21) for the case of 
unrestricted pion phase space, as defined by Eq. (63) 

where 

s na^a F + L F 1R + F 2R ^ ^ = ^ _ ^ ^ ^ = ^ _ ^ 

Because the factor a(q 2 )/AE 2 enters into the last term at the right-hand side of Eq. (64) 
also, the total RC in this case has a factorized form as well. 

There exists one more contribution caused by two hard photon emission when neither photon 
is emitted within the narrow cone along the the electron beam direction but, nevertheless, the 
collinear condition (4) is satisfied. This contribution cannot be calculated by the QER approach 
and has to be evaluated by other methods. In particular, the double bremsstrahlung lepton 
current tensor can be taken in limit m — > 0. To our understanding, due to the strong constraint 
(3) on the event selection, the corresponding contribution is small enough and does not affect 
the IES cross section on the one per cent level. Nevertheless, the theoretical evaluation of this 
contribution should be done and we hope to compute it elsewhere. 



6 Pair production contribution into the IES cross section 



The above considerations for the photonic radiative corrections to the IES cross section 
are appropriate if e + e~7r + 7r~ final states are excluded from the analysis. If not, there is an 
additional contribution caused by hard initial-state radiation with e + e~ pair production ]16] . 
The main part of this contribution arises due to collinear kinematics. In the framework of the 
NLO approximation, where only logarithmically enhanced terms are kept, the corresponding 
cross section can be written as 

, e+e~(c) / o\ o r e+e _ (c) , e+e _ (c) 
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where the functions P\{z) and ^2(2) can be extracted from the corresponding results for small- 
angle Bhabha scattering cross-section, given in Ref. [5^]. We present them here for complete- 



ness, i.e. 



P*(z) 



107 136 



1 + z 2 (1- z)(A + lz + Az 2 



3(1 -z, 



9 



— z z 

9 3 



+ 



\zT + 5z-7 
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In z 



+ 

A 

3z 

2 



62 



+ (1 + z) In z 



— r + - -4z 2 + hz 

9(1-2) 3 V 

In 2 2 + 4(1 + 2) lnzln(l - 



1 + 



z(l-z 
2(32 2 - 1 



ln(l - 2) 
■Li 2 (l - 2 



1-2 'V'/ v / x _ z 

Within the NLO accuracy, one has to compute also the contribution caused by the semi- 
collinear kinematics of the e + e~ pair production, when the final state electron belongs to the 
narrow cone along the electron beam direction while the positron does not. The corresponding 
part of the leptonic tensor was derived in [jST], and has the form 



r e+e- 



a 



d 3 k + dx(l + x 2 ) 



Air 3 £ + (l — x) 2 svi 



■{[ 



s 2 + vf 



+ 



q 2 v 2 



x 



g^u + 



2q 2 



1 



x 



where k + (e + ) is the 4 - momentum (energy) of the non-collinear positron, x is the energy 
fraction of the collinear electron and v\% = — 2p li2 A+. 

For unrestricted pion phase space, the differential cross section can be written as follows 



da e F +e - {s) 



AE 2 \2irJ 



a \ 2 dVLdx (1 + x 2 )[(l 
— I Lq- 



X 



z) 2 + z 2 ] 



X) 



(67) 



Neglecting terms of order m/E, we can integrate over the region shown in Fig. 2 with the 
substitution u/E — > x. In addition to the list of integrals (48), it is necessary to compute the 
following ones 



dx dQ 



zfi 1 
+ 



[l-xfVL z 2z(\ + z£) 2z 



1 + ln 



l + zi 



1 + z 2 
2z 2 



In z 



dx dQ 



:i-x) 4 a 



[l-z)(3 + z) 



3z 3 

l + z 3 
3z 3 



+ 



:i-z) 



(l-z) 2 (2 + z) 

Qz 3 (l + z£) 2 3z 3 (l + z£) 



(68) 



1-2 3 , l + z£ 
m 2 in ■ 



3z 3 



Using integrals (48), (68) and the definition (67), the contribution of the semicollinear 
kinematics of pair production into the IES cross section can be written in the form 



da 



e+e— (s) 



a 



dq 2 



a \ 



AE 2 \2irJ 



LoS(z,0 



(69) 



:i-z) In 



e 2r , {l-zfi , 2{l + z + z 2 )^ 1+zC 



1 + e 



1 + 



;i + zty- 



+ 



In 



2(1 + 2) 



1-2, 



ln£- — (1 -2 + 2 2 ) )hxz + Li 2 {-z£) + Li 2 ( ) -Li 2 (l - z) - Li 2 {-£) 
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Figure 4: The radiative correction to the IES Born cross section caused by e + e -pair production 



The contribution of e + e pair production into the IES cross section for the case of unre- 
stricted pion phase space now reads 



da\ 



dq 2 



a(q 2 ) a 
AE 2 2^ 



P(z,L )8 e p e , 8 e p e = 



a P 1 (z)L 2 + (P 2 (z) + S(z,O)L 



2tt 



P(z,L ) 



(70) 



with the function <5^ +e shown in Figure 4. In the case of restricted pion phase space, the 
differential distribution over the pion squared invariant mass can be written, in analogy with 
Eq. (55), in the following form 



da 



e+e (s) 
R 



dq 2 



12a(g 2 )/«\2 (1 + x 2 )dxdttT p Hp d{^ + - 

~Je 2 ~^\2^j (i- x yn z 2~ 27 



(71) 



where the integration region in x and Q is the same as in (67). Here, we use the following 
notation 



X2 



T 



2 \ Z +0--X-Z) ]+XX2 + (l~X - Z)X+ + X2 + X+ , X2 



X2 



2z{AK - Bf 



AE 2 ' 

A, B = xc__ — (1 — x — z)c+- , 



X+ 



2k + p. 
AE 2 



K(A 2 - B 2 ) 2 

c + _ = c + c_ + s + S- cos (<p + — if J) , c + = cos^ + , s + = sin9 + , 

and 6 + , <p + are polar and azimuthal angles of the final non-colinear positron. 

For the further integration on the right-hand side of Eq. (71) with respect to the angular 
pion phase space, one must determine the upper limit of c_ variation. It depends on q 2 , x, Q 
and cos (<p + — </?_) and can be obtained as a solution of the equation 



-|p+|c m + |p_|c_ + E(x + (1 - x - z)c+) = , 



(72) 
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taking into account that 



2Ez(A-KB) 

E_ = ^ = '- , E + = E(l +z)-E- . 

A 2 -B 2 

The total contribution of e + e~-pair production into IES cross section for the case of re- 
stricted pion phase space 

da e + e ~ _ da e + e ~ {c) da e + e ~ {s) 



+ (73) 



dq 2 dq 2 dq 

has no singularity and can be calculated numerically. 

We have considered in the above the contribution of kinematic regions where at least one 
collinear photon or an electron-positron pair is radiated by the initial electron. DA$NE con- 
ditions allow to select and detect also the same events when the collinear particles are emitted 
by the initial positron. Therefore, all the cross sections derived above have to be doubled. 



7 Conclusion 



The success of precision studies of the hadronic cross section in electron-positron annihi- 
lation through the measurement of radiative events [jn|, ^1], ^] relies on the matching level of 
reliability of the theoretical expectations. The principal problem is the analysis of radiative 
corrections corresponding to realistic conditions for event selection. 

In previous work |50] we discussed briefly an inclusive approach to the measurement of the 
hadronic cross section at DA$NE in the region below 1 GeV by the radiative return method, 
for the case in which the radiated photon remains untagged. This approach requires an exact 
knowledge of the final hadronic state and a precise determination of its invariant mass. Some 
additional constraints have to be imposed to make the detection of the ISR photon redundant 
and to avoid any uncertanties in the interpretation of the selected events. These additional 
constraints imply also the precise measurement of the pion 3-momenta. The KLOE detector 
at DA$NE offers a very promising possibility to realize such an inclusive approach to the 
scanning of the hadronic cross section by the ISR events. 

In this paper we compute the corresponding ISR Born cross section and the radiative correc- 
tions to it in the framework of the QRE approximation. It is shown that this approximation is 
quite appropriate, even at the Born level, and provides high accuracy for the IES cross section. 
The cases of unrestricted and restricted pion phase space are considered. In the first case the 
photonic contribution to the RC is calculated analytically with the NNLO accuracy and the 
contribution caused by e + e~-pair production within the NLO. The photonic RC is large and 
negative in a wide range of pion invariant masses. The physical reason for such behaviour of 
the photonic RC is very transparent: the phase space of additional real photon is restricted 
considerably by the constraints (3) and (4), and the respective positive contribution cannot 
compensate the negative contribution due to the virtual correction. The large absolute value 
of the first order photonic correction indicates unambiguously that the second order RC has to 
be evaluated. Moreover, the increase of the soft part of the RC (when z grows and approaches 
unity) requires summation of the leading RC to all orders. 

If the entire phase space for photons and e + e~-pairs is allowed, this problem is solved by the 
ordinary Drell-Yan-like representation in electrodynamics |32], |33|] with the exponential form 



22 



of the electron structure functions |28], |34| . For the tagged collinear photon events without any 
constraints on the phase space of additional particles, the corresponding representation was 
derived in Ref. |L6|], but the case considered here requires special investigation because of two 



non-trivial constraints (3) and (4) on the event selection. We hope to consider this problem 
elsewhere. 

The RC caused by e + e _ -pair production is positive and small, as compared with the absolute 
value of the photonic correction. Only in the region near threshold (small z), where the cross 
section is very small, it approaches approximately the same value. So, we conclude that the RC 
due to pair production described by Eq. (70) is adequate, and it must be taken into account to 
guarantee the one per cent accuracy. 

In the case of restricted pion phase space, we derived the analytical form of the acceptance 
factor for the Born cross section and the part of RC which includes contributions due to vitrual 
and real collinear photons and e + e~ pair. Even the choice of the detected pion angles between 
20° and 160° gives a value of the acceptance factor very close to unity, providing good statistics 
at all values of the hadron invariant mass. Concerning the contribution into the IES cross section 
caused by the semicollinear kinematics for the double photon emission and pair production, the 
respective acceptance factor cannot be calculated analytically. In this case formulae suitable 
for numerical calculation are given. 

One of the advantages of the IES approach discussed in this paper is a considerable decrease 
of the FSR background. The corresponding contribution into the IES cross section is suppressed 
by factors of order 9q due to the collinear restriction (4) on event selection, therefore we can 
ignore any RC to the FSR events and evaluate this background only at the Born level. The 
same is valid for the contribution caused by the ISR-FSR interference. Note that, for the tagged 
photon setup, the corresponding background is quite large and, in order to obtain the one per 
cent accuracy |Tj|, one needs evaluate the radiative corrections to it. If events with e + e~7r + 7r~ 
final state cannot be excluded by the experimental selection, the background due to the double 
photon mechanism of 7r + 7r~-pair production has to be evaluated as well. To our understanding, 
this background contains the same suppression factor 9q. In addition, restriction (3) selects 
very specific kinematics and does not allow to reach the region where the virtualities of both 
intermediate photons are small, and the corresponding cross section is the largest. In fact, due 
to this restriction, at least, one of the photons becomes far off-shell (with virtuality of the order 
—q 2 ). Thus, we expect that double photon mechanism contributes at the level of the RC to 
the FRS and cannot affect the IES cross section at the one per cent level. 
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Trust for a Fellowship. N.P.M. thanks INFN and Parma University for the hospitality. G. P., 
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